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Abstract 

We have investigated the behavior of three curvature invariants for 
Schwarzschild, Reissner-Nordstrpm, Kerr, and Kerr-Newman black holes. We 
have also studied these invariants for a Schwarzschild-de Sitter space-time, 
the 7 metric, and for a 2-\-l charged dimensional black hole. The invariants 
are h = , h = and I 3 = For 

all but the Kerr-Newman case these invariants serve as either horizon or sta¬ 
tionary surface detectors. The Kerr-Newman case is more complicated. We 
show that Ii vanishs on the horizon in any space-time with a Schwarzschild 
like metric. 

I. INTRODUCTION 

An outstanding challenge in the calculation of gravitational collapse to black hole con- 
hgurations is the identihcation of emerging black holes in the data on space-like slices of 
the space time Since the dehnition of event horizons involves knowledge of the future 

evolution of the data, the computational situation is not as clear-cut as one would like. In 
fact, the computational burden is shifted towards the identification of apparent horizons, 
as discussed in the previously cited references. Hoping to alleviate some of this burden, 
we have tried to develop criteria, involving the curvature tensor and its derivatives, which 
would indicate the existence of black holes. These criteria would, ideally, be complementary 
to those involving marginally outer-trapped surfaces. These criteria also would indicate a 
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trend in the data towards the formation of black holes. This hope has been fulhlled only 
partially. We take the opportunity to correct a result which has previously appeared in the 
literature. 

It has long been believed that an observer falling toward a black hole cannot tell when 
he or she has crossed the event horizon. This belief is based on the fact that the horizon is 
a global property of the geometry. The invariant 

h = ( 1 . 1 ) 

involves a covariant derivative of the curvature tensor and thus may be used to probe a neigh¬ 
borhood of the event horizon. This means that it has a chance of functioning as a horizon 
detector. This fact was hrst noticed by Karlhede, Lindstrom, and Aman who calculated 
Ji for the Schwarzschild, Reissner-Nordstrpm, Taub-NUT, and Kerr solutions. Karlhede, 
Lindstrom and Aman found that R vanished on the event horizon for the Schwarzschild, 
Reissner-Nordstrpm, and Taub-NUT solutions, and that it vanished on the stationary sur¬ 
face for the Kerr solution. Our results differ from those of Karlhede at el by a sign and due 
to our sign convention which is that of Misner, Thorne, and Wheeler [^. R has also been 
investigated by Tammelo and Kask [Q. To our knowledge, the invariants 

R = (1.2) 

and 

R = (1.3) 

have not been previously investigated. The curvature invariants are, in principle, locally 
measurable quantities [|,^|^, and thus could allow the local detection of an event horizon 
or a stationary surface. 

Using Mathematica and MathTensor, we have calculated R, R, and R for the Kerr- 
Newman geometry and for the special cases of Kerr, Reissner-Nordstrpm, and Schwarzschild 
black holes. We have also computed R for a Schwarzschild-de Sitter space-time, for the 7 
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metric IQ], and for a charged 2+1 black hole space-time. On a 300 MegHz PowerMac G3 
the calculations take several hours of run time for the longest (Kerr-Newman) cases. For 
vacuum space-times, Ji = I 3 and, of course, I 2 = 0. We conjecture that Ji serves as an 
event horizon detector for any static space-time where the horizon is regular. If the horizon 
is singular, it is not clear how to distingish the horizon from any other type of singularity. 

In addition to being intrinsically interesting horizon detection is of considerable interest 
in numerical relativity. The ablility to track the causal structure of a space-time is believed 
to be a crucial step in numerically evolving black hole space-times [|l|. For this reason 
considerable effort has gone into constructing apparent horizon hnders in numerical relativity 
The invariants discussed here maybe useful in that regard. 

II. STATIC AND STATIONARY SPACE-TIMES 

We have investigated the behaviour of Ji for both four dimensional space-times and a 
2+1 black hole space-time. 


A. Four Dimensional Space-Times 


1. Schwarzschild black holes 


We begin by discussing the Schwarzschild solution since it is the simplest case. 
For the Schwarzschild solution, the metric is 


ds^ = —(1 — 2M/r)dr + (1 — 2M/r) ^dr^ + r^{d 6 ^ + sin^ 9d(j) ), 


( 2 . 1 ) 


and we hnd that 


h = 


720 m2 (2 M - r) 


( 2 , 2 ) 


It is clear from equation (2.2) that the invariant is positive outside the event horizon, 
goes to zero on the event horizon, and is negative inside the event horizon. The invariant 
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has a maximum at r = 9M/4, and the only zeros of the invariant are at inhnity and on the 
event horizon. An in-falling observer can thus detect the presence of the event horizon by 
monitoring Ji. The in-falling observer may, for instance, measure the variation in oscillation 
frequencies from the equation of geodesic deviation 0. An in-falling observer can not save 
him or her self by turning back when the invariant vanishes for by that time it is too late! 
However, the invariant does provide a warning and the observer could turn back when the 
maximum of the invariant is reached. 


2. Reissner-Nordstr0m black holes 


The metric for the Reissner-Nordstrpm solution is 

ds^ = —(1 — 2M/r + Q‘^/r‘^)dt^ + (1 ~ 2M/r + /r^)~^dr‘^ + r‘^{d9‘^ + sin^ 9d(p^). (2.3) 

The Reissner-Nordstrpm case is somewhat more complicated, but the results for R are 
essentially the same as those for the Schwarzschild solution: R vanishes only at r = cx), on 
the horizon r_|_, and on the null surface r_. There is a positive extremal value of R just 
outside of r+. The full invariant is 

16 {7QQ^-108 MQ^r+ A5M^r^) {Q^ + r (-2 M + r)) 


R = 


U2 


(2.4) 


and we see that R = 0 at r = cx). The vanishing of R can be reduced to hnding the roots 
where the numerator of equation 2.4 vanishes. 

There are two hnite complex roots and two hnite real roots. The two real roots are the 
null surface r_ = M — y/M'^ — Q'^ and the event horizon r+ = M -|- — Q‘^. This result 

was hrst obtained by Karlhede, Lindstrom and Aman [Q, although there is a typographical 
error in their result for R for the Reissner-Nordstrpm case. Figure one shows a typical plot of 
R as a function of r in the vicinity of r+ and hgure two shows a typical plot of R as a function 


of r in the vicinity of r_. We notice that R is positive outside r+, negative in the region 
between r+, and r_ and 

, 2(9Q2+i\/l4Q2) 

and r = —- tvtt- -■ 


between r+, and r_ and positive inside r_. The complex roots occur at r = ^ ^ 


15 AT 
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For the Reissner-Nordstr 0 m case, I 2 

^ + ^ (- 2 M + r)) 

^2 -- -^2 -• 

We observe that, like Ji, I 2 vanishes at r+ and r_, bnt that nnlike Ji, I 2 has no complex 
roots. Figure three shows I 2 in the region of the horizon. 

The invariant is given by 

48 (Q2_ 2 Mr+ r2) (2236MgV +15 m 2r2) ^ ^ 

^3 =- ^2 -• ( 2 - 6 ) 

Like 1 1 and / 2 , Is vanishes at r+ and r_. In addition J 3 has two complex roots which he at 
^ and r = The complex roots of h and Is are not the same but 

near the horizon R and are numerically close in the sense that |(/i —/ 3 )//i| 1 . Figure 
four shows the difference between Ji and in the vicinity of the horizon for a particular 
choice of parameters. 

For an extremal black hole Q = M the two roots at r+ and r_ coalesce into a double 
root at r = M. If > M^, then Ji, / 2 , and J 3 have no real roots. 


3. Schwarzschild-de Sitter spacetimes 


For the Schwarzschild-de Sitter space-time the metric is 

ds^ = —(1 — 2M/r -I- Ar‘^/3)dt‘^ + (1 ~ 2M/r -|- Ar'^/3)~^dr'^ + r‘^{d6‘^ + sin^ 9d(j)^), (2.7) 


and 


h 


240m2 (- 6 M + 3r-r3A) 


( 2 . 8 ) 


This invariant is zero when 1 — 2M/r — Ar‘^/3 = 0, which is the event horizon. Thus, Ji 
seems to serve as a cosmological event horizon detector. 
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Space-Times with Schwarzschild like metrics 


We call a Schwarzschild like metric one with the line element 

ds"^ = —A(r)dt^ + A{r)~^dr'^ + r'^dO'^ + sin^ 9d(()^. (2.9) 

where A{r) is an arbitrary function of r. For this metric we can calculate Ji for a general 
A{r) and we hnd that 

Ji = A{r) [ (32 + 32 A{rf + 16 A!{r-f - 32 A{r-) (2 + r W(r)) 

+8 W'(r)^ — 16 r W(r) 2 + W'(r)^ + r® /r®, (2-10) 

where the primes indicate differentiation with respect to r. We notice that Ji vanishes when 
A{r), the norm of the time-like Killing vector, vanishes. In addition Ji will vanish when 

32 + 32 A{r-f + 16 A\rf - 32 A{r) (2 + r W(r)) 

+ 8rM"(r)^ - 16rW(r) (-2 + rM"(r)) + r® W"(r)^ = 0. (2.11) 

With the appropriate choices for A{r), Ii reduces to either the Schwarzschild or the Reissner- 
Nordstrpm result. We have also calculated the invariants I 2 and for this space-time and 
hnd that 

I 2 = A{r) [ 2 A -\- 24 A{r)^ + W(r)^ -|- 24 A"{r) + 10 W'(r)^— 

2AA{r-) {2 + A!'{r)^ + Ar^ A!'{r) A^^\r) + X” (r)^— 

4rM'(r) (2kl"(r)+ rkl"'(r)))]/2r® (2.12) 

and 

J 3 = A{r) (^40 -|- 40 A(r-)^ -\- 40W(r)^ — 40 X\r) -\- 10W'(r)^ -|- 8 A'"(r) — Ar^ A''(r) A"'{r)-{- 
X” ij-)^ + S r X{r) (^10 — 5 A"(r)-f-R"'(r)^ — 8 R(r) (lO-l-lOrW(r) — 

5 X'(r)X” /3r®. (2.13) 

If the space-time is a vacuum space-time then the held equations imply 
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1 — A{r) 


(2.14) 


A'(r) = 


r 


and I 2 = 0, while 


180 (-1 H- /l(r))^yl(r) 


(2.15) 


5. Kerr Black Holes 

The metric for the Kerr solution is 

ds^ = + d9‘^ + (r^ + a^) sin^ Odcfy^ — dt^ + (^asin^ 9d(j) — dt^ , (2.16) 

where p‘^{r,9) = + a^cos^^ and A(r) = — 2Mr + a^. For Kerr black holes the only 

invariant of interest is Ji, and we find the remarkably simple result that 

Ji = 368640 (^—2Mr + cos(6*)^^ x 

— 28 r® cos(6*)^ + 70 cos(6*)'^ — 28 a® cos(6*)® + a® cos(6*)® 

(a^ + 2 cos(2 0))® 

The form of our result is different from the one given by Karlhede at el [Q but the two agree 
up to the previously mentioned sign convention. 

For the Kerr solution, the invariant Ji has 10 roots. Two of these involve the mass of the 
hole and are at the inner and outer stationary limits given by = M ± cos^ 9. 

The other eight roots depend only on a. The fact that Ji vanishes at the stationary limit 
rather than the horizon is not surprising since the norm of the Killing vector ^ vanishes on 
the stationary surface. The Killing vector ^ ,where D.h = ^2 is tangent to the 

null geodesic generators of the horizon, p. 

The eight Kerr roots that depend only on a may be expressed in a simple form and are 
given in Table 1. We observe that the roots are paired; for every root at a positive value of 
r there is a root at the corresponding negative value of r. Each root in the left column of 
Table 1 is paired with the corresponding root in the right column. The existence of roots for 


(2.17) 
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negative values of r is perhaps not surprising since the maximally extended Kerr solution 
includes both positive and negative values of r H- As a —> 0, the mass independent roots 
coalesce to zero. If then the two mass dependent roots become complex. 

The zeros of Ji are shown in hgure 5. The behavior of Ji in the vicinity of the horizon is 
shown in hgure 6. 


6. Kerr-Newman Black Holes 


The most complicated case we consider is the Kerr-Newman black hole. The metric for 
the Kerr-Newman solution is 


ds^ =—{/S./ p^) {dt — 6d(j^ + 6/p^) dcj) — adt 

+ {p^//\)dr^ + p‘^d9‘^, 


(2.18) 


where p^(r, 9) = r'^ + cos^ 9 and A(r) = — 2Mr + a'^ + Q'^. 

For the Kerr-Newman solution, all three invariants are distinct and in general do not 
vanish on either the horizon or on the stationary surface. We hnd that 


h = 


8192 (r® (76 Q® - 260 M r -h 26176 - 90 

-lOSMQ^r^+ 45M2r^(44Q^ - 36MQV)) + x 

(-744 Q® -h 3540 M r - 5364 ^2 _ 2520 1980 MQ^r^- 

1215 (-156 -h 324 M Q'^ cos{ef + 6a‘^r'^ (98 Q® - 706 M r -h 1545 ^2 

-1050 MW^- 222 315 (-26 -h 30 M r)) cos{ef+ 

2 a® (-16 Q® -h 302 M r - 1170 ^2 g4 ^2 ^260 - 1350 M + 945 r*+ 

2 ad- (llQ^-45 cos(0)® - o® (76 + 90 r - 720 M gV - 45 (g2-27r2)) cos(0)^+ 

45a^°M2cos(0)^°)] /(a2 + 2r2 + a2 cos(20))®, (2.19) 


/2 = 


—2048g"^ (— (r^ (4a^ -|- 5 (Q^ — 2Mr-\- r^))) — (4a^ — -|- 2Mr) cos(0)^ 4-5 0 “^ cos(0)"^^ 

(a^ -|- 2 r^ -|- cos(2 0))”^ 


( 2 . 20 ) 








and 


h = 


24576 (r® (22 Q® - 80 M Q V + 87 ^2 22 - 30 - 36 M + 15 /+ 

12 (g^ - M g V) ) + a (-254 g® + 1192 M g V - 1788 g^ - 244 g^ r2+ 

840 + 660 M - 405 (-60 g^ + 108 M r'^'^ cos{ef+ 

2 (97 g® - 704 M gV + 1545 g^ - 1050 - 222 M g^ + 315 r^- 

30 (g^ -UQ'^r^j cos{ef + 2 a® (-5 g® + 100 M gV - 390 + 122 g^ r^+ 

420 - 450 M Q"^ + 315 M + 6 (q‘^ -bM g^ r) ) cos(6»)® - a® (22 g^+ 

30 mV - 240 M gV - 15 M^ (g2 - 27 ) cos(6»)V 

15a^°MVos(0)^®))] /(a2 + 2r2 + a2 cos(20))^ (2.21) 


Although Ji and /a are distinct, the surfaces on which they vanish have, for typical values of 
the hole parameters, similar shapes. Both of these invariants have 10 roots, of which two depend 
on the mass of the hole. The other eight roots depend only on the charge and angular momentum. 
The surfaces on which Ji vanishes are shown in figures 7, 8, 9, and 10. The roots of Is are similar. 
As a ^ 0, the roots move toward the Reissner-Nordstrpm case. The physical significance of the 
mass independent roots is unclear. As a ^ 0, the two mass dependent roots move to the inner 
and outer horizons and the mass independent roots coalesce into the two complex roots of the 
Reissner-Nordstrpm case. This is shown in figure 11. 

The invariant I 2 has four roots which are shown in figure 12 for a hole with M = 2,g = 1, 
and a = 1. We notice that the outer root lies in the ergosphere, and that the second root weaves 
in and out of the inner horizon. The two remaining roots lie inside the inner horizon. In the limit 
a ^ 0, the two inner roots coalesce to zero. The roots of I 2 for an extreme Kerr-Newman hole 
are shown in figure 13. Although it is not obvious from the figure, I 2 also vanishes on the event 
horizon at 0 = 0 and vr. The two mass dependent roots remain real for some values of 9 even for 
the ultra-extreme case where > M^. In general none of the invariants vanishes on the 

horizon nor on the stationary surface, except at 0 = 0. Some of the roots of Ii, I 2 , and I 3 are 
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complex for some values of a, Q and 9, as they must be, since both Ii and I 2 have two imaginary 
roots for the Reissner-Nordstrpm case. 


7. The 7 metric 


A static, axially symmetric, asymptotically flat vacuum space-time is provided by the 7 metric 
1^], a two parameter Weyl solution. The 7 metric serves as an example of a space-time where the 
horizon, depending on the values of the parameters M and 7 , may be regular or not. The 7 metric 


IS 


,2 /, 2M\^,2 a 2 M \-^ 

ds^ = -{1 - 1 dt^ + [1 - 


— 2Mr 


7-1 


— 2Mr) 


— 2Mr + M sin 
—dO'^ + (r^ — 2Mr) sin O’^dcjP' 


dr^+ 


( 2 . 22 ) 


(r^ — 2Mr + sin 6 ‘^y 

and corresponds to a total mass of 7 M. The 7 metric reduces to the Schwarzschild metric when 
7 = 1. While the complete analysis of the 7 metric is rather complicated, we mention only a few 
of the features which have a bearing on the present calculations. For 0 < 7 < 2 , 77^1 the null 
surface r = 2M is singular along the symmetry axis 0 = 0 or vr. For 7 > 2 the null surface r = 2M 
is regular along the symmetry axis. This can be seen by looking at Iq = . We have 

computed both Iq and R for general 9 but the results are long. Along the symmetry axis we find 

(2.23) 


Iq = 


48M2(i-^) ^-f^{M-r + M-fY 
r*{-2M + r)^ 


and 


80 m 2 (3r2-6Mr {l + 'y) + 2M^ {2 + 3 j +'y‘^)y 


h = 


(2.24) 


r^{-2M + rf 

which reduces to the Schwarzschild result if 7 = 1 . R vanishes on the horizon if 7 > 2 . There 
are also two zeros of R found symmetrically about r = M(1 -|- 7 ). In short R, is singular on the 
horizon when the horizon is singular and vanishes on the horizon when the horizon is regular. The 
belief that Weyl solutions are important in collapse situations hinges on the hope that reasonable 
interior solutions can be found for them. 
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B. 2+1 Black Holes 


We have also computed Ii, I 2 , and Is for the BTZ 2+1 dimensional black hole. The charged 
BTZ black hole found by Bahados, Teitelboim, and Zanelli 1]^ has a metric given by 


-1 


ds^ = — \ ^ — M + Inr \ dt‘^ + I ^ — M + Q"' In r 1 + r'^dcj)^. 


cr^ 




(2.25) 


For the charged BTZ black hole 


12 + + 4 cj^ + 3 Q‘^ 

£j4 


/n = 


(2.26) 


The scalar curvature is 


R = 




(2.27) 


and the invariant Ji is 


^ 20 (_M + (r/u)2 + Q2 iog(^)) 

J.6 ’ 

which is seen by inspection to vanish on the horizon. The invariant I 2 is found to be 

6 (r^ — M + g^ cr^ log(r)) 


(2.28) 


(2.29) 


The space-time is Weyl flat, so I 3 = 0. 

Since the geometrical and physical interpretation of the curvature invariants is unclear it may 
be useful to rewrite I 2 in terms of Using the Einstein equations, we find that I 2 may be 

written in terms of invariants involving and R. It is straightfoward to show that 

I 2 = 4: (^emi — + -^R-\R'^, (2.30) 

where 


IeM 2 = . 


(2.31) 


(2.32) 
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We find that 


I EMI = 


2 (r^ — M a'^ + Q"^ a'^ log(r)) 


(2.33) 


IeM2 = 


4 (r^ — M cT^ + log(r)) 


r® (j2 


(2.34) 


and 


^ 4Q^ (r^ - Mcr^ + Q2cr2 log(r)) 

n.\W — - 


r® cj2 


(2.35) 


Both the electromagnetic invariants and the scalar curvature invariant vanish on the horizon, and 
are in fact proportional to l 2 - 

The uncharged black hole space-time has constant scalar curvature and is, for our proposes, 
uninteresting since Ii is identically zero. 


III. CONCLUSIONS 

The invariants Ii, I 2 , and Is, which are locally measurable, serve as horizon or stationary 
surface detectors for a variety of space-times but not for Kerr-Newman black holes. We have 
shown that Ii serves as a horizon detector for any Schwarzschild like space-time and conjecture 
that Ii serves as a horizon detector for regular horizons in any static axially symmetric space-time. 

It is not known if an invariant can be found that will vanish on the horizon in the Kerr-Newman 
case. 

A problem of considerable current interest in numerical general relativity is the detection of 
emerging horizons in collapse problems. The invariants discussed here, particularly Ii and I 2 , may 
be useful tools for the numerical detection of such emerging horizons. 
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TABLES 



r = ^—1 + a cos{6) + \JA — 2 \pi \J a? cos(0)^ r = — « cos(0)^ — — 2 cos(0)^ 

r = — ((—1 + \/2^ a cos(0)'j + ^4 — 2 \/2 cos(0)^ r = ^—1 + \/2) a cos{6) — \JA — 2 \f2 \Ja? cos(0)^ 

TABLE 1. The eight roots of Ji = 0 that are independent of the mass of the hole for the Kerr 


solution. 
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FIGURES 



FIG. 1. Ii for the Reissner-Nordstr0m solution as a function of r for the case M = 2, Q = 1 in 
the vicinity of r+ 
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FIG. 2. Ii for the Reissner-Nordstr0m solution as a function of r for the case M = 2, Q = 1 in 
the vicinity of r_ 



FIG. 3. I 2 for the Reissner-Nordstr0m solution as a function of r for the case M = 2, Q = 1 in 
the vicinity of r_|_ 


16 




FIG. 4. Is — Ji for the Reissner-Nordstr0m solution as a function of r for the case M = 2,Q = 1 
in the vicinity of r_|_ 



FIG. 5. The roots of Ii for the Kerr metric. The inner and outer stationary surfaces are 
shown in black. There is a root on each stationary surface. Four other roots are shown in grey. 
The remaining four roots are shown as dashed lines with a different dashing used for each root. 
Two of the roots are inside the inner stationary surface and are thus hard to see. All the roots 
except the two which lie on the stationary surfaces depend only on a. 
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Il 



FIG. 6. Il or the Kerr metric as a function of r for the case M = 2,a = 1. The approach is 
along the 0 = 0 axis and the stationary limit is indicated by the black dot. 



rcos "d) 


FIG. 7. The outer roots of Ii for the Kerr-Newman metric with M = 2,a = 1,Q = 1. The 
horizon and outer stationary surfaces are shown in black. The three outer roots are shown in gray 
and as dashed lines. 
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r sin(d) 


M 



FIG. 10. The other three inner roots of Ji for the Kerr-Newman metric with 

= 2, a = 1, Q = 1. These roots all lie outside the inner stationary surface. 
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a = 1 


a = 1 / 2 


r sinC-d) r sinC-d) 



a = 1 / 4 


a = 1/ 16 


r sin(i3-) r sin(i3-) 



FIG. 11. The three outer roots of Ii for the Kerr-Newman metric as a ^ 0 with M = 2,Q = 1. 
As the hole spins down the roots move toward the Reissner-Nordstrpm case. The outer horizon 
and stationary surface are shown in black. The two mass dependent roots move to the inner 
and outer horizons and the mass independent roots coalesce into the two complex roots of the 
Reissner-Nordstrpm case. For reasons of scale only one of the mass independent roots is shown. 
The root shown with the dashed line becomes complex as a ^ 0. 
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r sin( 'O') 



cos( 0) 


FIG. 13. The t'wo outer roots of I 2 for the Kerr-Ne’wman metric for an extreme case ■with 
M = 2, a = 2^/^, Q = 2^/^. The two outer roots have joined each other. The light gray line shows 
the same root as the outer gray line in figure 12. The dashed line is same root as the inner gray 
root in figure 12. The horizon and the stationary surface are shown in black. 
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